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Abstract 

The most general exclusion single species reaction-diffusion models with 
nearest-neighbor interactions one a one dimensional lattice are investi- 
gated, for which the evolution of full intervals are closed. Using a gener- 
ating function method, the probability that n consecutive sites be full is 
investigated. The stationary values of these probabilities, as well as the 
spectrum of the time translation generator (Hamiltonian) of these are ob- 
tained. It is shown that depending on the reaction rates the model could 
exhibit a dynamical phase transition. 
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1 Introduction 



Different methods have been used to study non equilibrium statistical models; 
including analytical and asymptotic methods, mean-field methods, and large- 
scale numerical methods. Although mean field solutions may be suitable for 
higher dimensions, in low dimensional cases fluctuation effects may cause mean 
field results to differ from the real ones. One dimensional models, which are 
usually easier to investigate, helps us gaining more knowledge on systems far 
from equilibrium [TMl3j. Many techniques are used to obtain exact results in 
one-dimensional models. The empty interval method (EIM) is one of them. 
The full interval method (FIM) is another, which is equivalent to EIM. In EIM, 
the aim is to calculate the probability that n consecutive sites be empty, E n . 
In FIM, one studies the probability that n consecutive sites be full, F n ,). In 
[14H17] , one dimensional models have been studied in which some of the reaction 
rates are infinite, so that it is the (finite) diffusion rates which determine the 
evolution. These models have been through EIM. In [18], a system with three- 
site interactions has been studied by EIM. In [19] . a generalization of EIM has 
been used to investigate a one dimensional Potts model with g-states, in the zero- 
temperature limit. In [20) . a classification has been given for one dimensional 
reaction-diffusion models with nearest neighbor interactions, which are solvable 
through EIM, by which it is meant that the evolution equations for E n 's are 
closed. In such systems, the time derivative of the empty intervals is linear in 
the empty intervals (and contains no other dynamical variable). If some extra 
conditions are met (there is no reaction which produces particles in two adjacent 
sites) the evolution equation for E n becomes n-independent. This makes solving 
the evolution equations easier. In [21H23) . these were generalized to systems with 
multi-species systems and multi-site interactions. In [23], the extra condition 
that the evolution equations for E n be n-independent was relaxed, and solvable 
(in the sense of EIM) models on continuum were studied. 

In [25], a generalization of EIM was introduced, by which a model was 
investigated which was not solvable through conventional EIM. In |26,27J, EIM 
has been used to study the coagulation-diffusion process on a one dimensional 
chain. 

In a recent article, the most general exclusion single species one dimensional 
reaction-diffusion models with nearest-neighbor interactions were studied, which 
are both autonomous and can be solved exactly through the full interval method 
[25] , There, using a generating function method, the general solution for F n , the 
probability that n consecutive sites be full, and some other correlation functions 
of number operators were explicitly obtained. 

In this article we relax the condition of autonomy, and study the most general 
exclusion single species reaction-diffusion models with nearest-neighbor interac- 
tions on a one dimensional lattice, which can be solved exactly through full 
interval method. The change of empty interval to full interval is, of course, not 
important, as a simple interchange of particles and holes would do that. The 
scheme of the paper is as follows. In Section 2, the most general exclusion single 
species one dimensional reaction-diffusion models with nearest-neighbor inter- 
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actions are introduced, which can be solved exactly through FIM. In Section 3, 
a generating function method is used to obtain the time independent solutions 
to the evolution of full intervals. In Section 4, a similar generating function 
method is used to study the spectrum of the time translation generator (Hamil- 
tonian) of the full intervals, it is shown that the model may exhibit a dynamical 
phase transition. 



2 Full interval equation 



Consider a one dimensional lattice, any site of which is either occupied by a 
single particle or empty. The full interval F n is defined as the probability that 
n consecutive sites be full 

n 

F n := P(T^~i). (1) 

An empty (occupied) site is denoted by by o (•). Assuming only nearest neighbor 
interactions, it can be seen ( [24]) that the necessary and sufficient condition 
that the evolution equations for F n 's be closed is that onl the following reactions 
be present, and their rates be related to each other as follows. 



oo — >• 




(2) 



The equality of some rates, means that the rates of the reactions (oo — > »o) and 
(o» — > »o) are the same, and the rates of the reactions (oo — » o») and (»o — > o») 
are the same as well. 

As in, for example |24) . it can be seen that the time evolution equation for 
F n is 



^ = (r 1 + r 2 )(F„_ 1 



F n+1 -2F n )- (<zi + q 2 ) (F n - F n+1 ) 
(n - 1) (wi + w 2 +w)F n - (wi + w 2 + 2w) F n +i, n>2 



-T7 1 = (ri + r 2 ) (1 + F 2 - 2*i - ( Ql + q 2 ) (F, - F 2 ) 
at 



— (u>i + w 2 + 2 w) F 2 . 
It is seen that defining 



F 



(3) 
(4) 
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the evolution equation for Fi takes a form similar to that of other F ra 's: 
dF 

= (ri + r 2 ) (F„_x + F n+1 - 2 F„) - (ft + g 2 ) {K - F n+1 ) 

— (n — 1) (»i + W2 + F„ — (wi + w 2 +2w) F n +i, n > 1. (5) 

Comparing this with the similar expression in, it is seen that in |24j there is no 
term analogues to the last term here. In 28 , it is assumed the models to be 
autonomous, which leads to removing the term F n+ i from (|5j). 
Defining 

t := (ioi + w 2 + w) t, (6) 



one arrives at 



where 



dF 

= bF n - 1 -(a + n-l)F n + cF n+1 , (7) 



2 (ri + r 2 ) + qi+q 2 

a := 



W + Wx + W 2 

r\ + r 2 



W + Wi + w 2 ' 



n+r 2 + qi+q 2 -2 w - w 1 - w 2 

c:= . (8) 

W + Wl + w 2 

The case (wi + w 2 + w) = 0, will be dealt with separately. For simplicity, 
hereafter the symbol t is used in place of t, so 

dF 

=6F„_x - (a + n-l)F n + cF n+1 . (9) 
A generating function F is also defined, which will be used later: 

n X )-.= f:^f. (10) 

n=0 

It is seen that 

F n =F^(0), (11) 
where F^ is the n-th derivative of F. 



3 The time independent equation 

The time independent solution F st satisfies 

b F°U - (a + n - 1) J* + c F r f +1 = 0, (12) 

which yields 

E £z Tv [ cF «Vi + bF n-i -(a + n-1) F?] = 0, (13) 

n=l \ >' 
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so that 



Defining 



d 2 F st dF st 



dx 2 



d.z: 



one arrives at 



So, 



F st (x) := z (1 ~ Q) G st (z), 
z := 2 y/b (c — x), 



z 2 + z ^_ + [z 2 _ (1 _ fl)2] gat = Q 



dz 2 



dz 



G s \z) = aJ Q _ 1 (z)+/3J 1 _ Q (z), 



(14) 

(15) 

(16) 
(17) 



where J„ is the Bessel function or order v, and a and /3 are constants. So, 

F st (x) = a [2^b(c-x)\^ J^y/bic-x)] 

+ (5 [2yjb (c - x)] (1 - a) h-a[2xjb (c - x)]. (18) 

As is in [0,1], the convergence radius of the series defing F st is infinity. 
So the the generating function F st as a function of x is analytic on the entire 
complex plain. Using 



3„(z) 



r(i/ + i) V2 



\z\ « 1, 



the analyticity of F st at (x — c) demands that 



As 

one has 



(3 = 0. 



r — L i 



F st {0) = 1, 



which can be exploited to obtain a. So, 

pst = [2y/b{c-x)]^ Z a -x[2^b{c-x)] 



(19) 

(20) 

(21) 
(22) 

(23) 



(2V&c)(i-«) J a _!(2V67) 
Using (fTTj) . one can find F n . One has 

1 / 26 d 

(2 VM^Ja-i^ v^) V~»~dl, 

(2 6)"z 1 - Q -"J a _ 1+ „(z), (24) 



_d_ 

da; 



F st (x) 



Z 1 — Ja-^Z)], 



1 



(2V6c)i-«J Q _ 1 (2v'6c) 
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which results in 

n/2 



F *t_(b\ n " J a -i +w (2V6c) 

" W J«_l(2V^) ' 

There are two limiting cases to be studied separately: 
i: b = 0. 

This can be studied as the limit b — > of the general case. Using (|25|) and the 
limiting behavior (Tl9|). one arrives at 

iff = J {a) , b n , b«l, (26) 



T(a + n) 



which leads to 



lim Ff = S^. (27) 

This is the case if a is nonzero. If a and both vanish, the result could be 
obtained directly from (fT2"|) to be 



(28) 




where 

0<P<^, (29) 

and it has been assumed that at least one of the rates are nonvanishing. One 
notes that p is in fact the density of the particles, and the restriction on its 
value results from the fact that in the stationary configuration no two-adjacent 
sites are full. 

These results are expected. If b vanishes but a does not, there are no reac- 
tions which produce particles but there are reactions which annihilate particles, 
whether the particles are adjacent to holes or other particles. So at large time 
the lattice becomes empty. If both a and b vanish, there are no reactions which 
produce particles, but there are reactions which annihilate particles, only if 
there are two neighboring particles. So particles will be annihilated, but when 
there are particles with empty neighboring sites, they will survive. 

In fact, regarding the stationary solution as the large-time solution, one can 
obtain an expression for p in terms of the initial conditions. As both a and b 
vanish, one has 

dF 

-^ = -{n-l)F n + cF n+1 . (30) 

Defining 

F n (t):=exp[(n-l)t]F n {t), (31) 

equation (1501) is recast to 

dt 



cexp(-t)T n+1 , (32) 
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which results is 

^w = E [c ' ce fc T ( " t)]fc ^ +fc (o), 

so that 

lfc 



kl 

k=0 



oo u 



(33) 
(34) 



F n (t) = exp[-(n - 1) t] £ [C CC ", P( t)] ' F n+fc (0) 

from which, 

fc=0 

or, 

OO 

fc=0 

One notes that when a and b both vanish, c is nonpositive. 
ii: wi + U>2 + w = 0. 

In this case all of the rates Wi , u> 2 , and w should vanish. One arrives at 

{ri + r 2 ) - (2 n + 2 r 2 + 9i + «&) + (n + r 2 + ft + g 2 ) F^ +1 = 0, 
F st = 1. (37) 

The solution to (|3"T)) is 

V + r 2 + gi + g2 / 

where 

< C < 1, (39) 
provided of gi and g 2 , at least one is nonvanishing. Otherwise 

F* = l, (40) 

which corresponds to a full lattice. 

4 Relaxation towards the time independent equa- 
tion 

To study the spectrum of the time translation generator (Hamiltonian) of the 
full intervals, again the generating function is used. Defining 

F^:=F n -F%, (41) 

one arrives at 

^ =bF*l l (a + n 1) + cF% 1 , (42) 
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with the boundary condition 

F dy = 0. (43) 

Equation (|42p is of the form 

dF dy 

-£- = (hF*) n , (44) 

where h is a linear operator. To find this relaxation time, one should obtain 
the eigenvalues of h. The eigenvalue with the largest real part, determines the 
relaxation time. Denoting the eigenvector of h corresponding to the eigenvalue 
E by ij)Et one nas 

ElpEn = bljjEn-l ~ (a + n - 1) IpEn + CtpEn+1, (45) 

where E is the corresponding eigenvalue. This is similar to (|12p . with a replaced 
by (a + E). So repeating similar arguments, one arrives at 

M*) = " (c - x)^- a - E V 2 J a+E -x[2Vb (c - x)}. (46) 

But now the boundary condition is 

MO) - 0, (47) 

which results in 

J a+E _!(2VM = 0. (48) 

This gives the spectrum of h. 

In the case b = 0, one can find more explicit forms for the eigenvalues and 
eigenvectors. Starting from (|45)) . one arrives at 



E + a + n - 1 
tpEn+i = WEn- (49) 

This recursive relations shows that ipEn tends to infinity as n tends to infinity, 
unless there is a k so that ipEn's vanish for n > k. This happens if E is equal 
to one of Ek's, where 

E k = l-a-k, (50) 

and k is a positive integer. Denoting the corresponding eigenvector by ipk instead 
of ipE, one would arrive at 

(-c) fc_ ™ 

^ = r(fe + i-n)^ fe - (51) 

The relaxation time of the system is obtained from the largest real part of the 
eigenvalues, which in this case is (—a). So 

t=~, (be) =0. (52) 
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Defining 



e := E + a - 1, 



(53) 



it is seen from (|4"8"|) that e depends on only the product (6 c) . So the expression 
(|5"P|) holds for the case c = as well. The case c = has already been discussed 
in greater detail in [25] 

From ([8]) is it seen that 6 is nonnegative, while c can change sign. These 
expressions also show that 

-l<(6c), (54) 

but there is no upper limit for (b c) . 

For (be) near zero, one can find the leading correction to (|5ip as follows. 
Defining 

5k ■= £fc + k, (55) 

(56) 

(57) 

(58) 
(59) 



equation (|4*8"f becomes 

J_ fe+St (2Vte) = 0, 
which is, up to the leading order, equivalent to 

(6c) fe 



1 



k\(k-l)\8 k 



showing that 



So one has, 



4 = 



(bc) k 
kl (k - 1)! 

1 



a — be 



0, 



|6c| < /c 2 



|6c| < 1. 



If (6 c) is positive, there is an inner product with respect to which h is 
Hermitian. So if (be) is positive, the spectrum of h is real. Increasing (be) from 
zero, the values of e^'s are also increased. For large values of (be) one can find 
the values of e^s using various asymptotic expressions of the Bcssel functions. 
One arrives at 



2 \/bc — <Xk(bc 



,1/6 



4V6< 



-2k 



2 

(bc) k 



k\(k-l)\ 



+ 



k < (6 c) 1 / 3 



|e fe |«(6c) 1 /2 



fc> (6c) 1 / 2 



(60) 



where (— ttfc) is the fc'th zero of the Airy function. The largest of e^'s determine 
the relaxation time. One then has 



1 

a- 1 + 2 Vbc 



(be) > 1. 



(61) 
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(be) 




Figure 1: 

The plot of (6 c) versus e, corresponding to the four largest eigenvalues of h 



From (jSJ is it seen that for nonnegative c, 

(a-l)>2Vbc. (62) 

For a negative (be), however, there are cases where the spectrum of h is not 
real. A plot of (6 c) in terms of (real) e for 

J e (2\Z&c)=0, (63) 

shows that there is one minimum for (6 c) for each interval e € (—2 n, —2 n + 1), 
where n is a positive integer, figure 1. However, only the minimum correspond- 
ing to e 6 (—2,-1) is larger than (—1). So there is a critical value for (6c), 
at which two of the eigenvalues of h (the largest and the next largest) become 
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equal, and if (b c) is less than that critical value, two of the eigenvalues of h 
become nonreal (complex conjugates of each other. Denoting that critical value 
of (be) by 7, and the corresponding value of e by e tr , one arrives at the following 
approximate expression for e and (be) near their critical value. 

(6c) = 7 + (64) 

where the numerical values of the constants in the above equation are 

7 = -0.401873, 
v = 0.754464, 

e c = -1.697524. (65) 

One then arrives at the following expression for e\, for (be) near the critical 

value 7: 

£1 = £ c + \]v(bc- 7), (66) 

resulting to 

( 1 

a - 1 - £ c - \Jv (be- 7) 

(67) 

, (bc)<j 

\, a - 1 - £ c 

So the derivative of the relaxation time with respect to (be), is infinite for 
(6c) — > 7 + , and zero for (be) — > 7". This model shows a dynamical phase 
transition. 
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